1 

Approaching the Capacity of Large-Scale 
MIMO Systems via Non-Binary LDPC Codes 

Puripong Suthisopapan, Student Member, IEEE, Kenta Kasai, Member, IEEE, 
Anupap Meesomboon, and Virasit Imtawil 

Abstract 

In this paper, the application of non-binary low-density parity-check (NBLDPC) codes to MIMO systems which 
employ hundreds of antennas at both the transmitter and the receiver has been proposed. Together with the well-known 
low-complexity MMSE detection, the moderate length NBLDPC codes can operate closer to the MIMO capacity, 
e.g., capacity-gap about 3.5 dB (the best known gap is more than 7 dB). To further reduce the complexity of MMSE 
detection, a novel soft output detection that can provide an excellent coded performance in low SNR region with 99% 
complexity reduction is also proposed. The asymptotic performance is analysed by using the Monte Carlo density 
evolution. It is found that the NBLDPC codes can operate within 1.6 dB from the MIMO capacity. Furthermore, the 
merit of using the NBLDPC codes in large MIMO systems with the presence of imperfect channel estimation and 
spatial fading correlation which are both the realistic scenarios for large MIMO systems is also pointed out. 

Index Terms 

Large MIMO systems, Hundreds of antennas, Non-binary LDPC codes, MMSE detection, MF detection, Soft- 
output detection. 

□DO 

I. Introduction 

The multi-input multi-output (MIMO) systems are the transmission systems that use multiple antennas at both 
sides of the communication ends (TJ, (2). In the past few years, more and more attentions have been devoted to 
large MIMO systems, e.g., MIMO systems with tens to hundreds of transmit/receive antennas Q, (4), Q since 
the large MIMO systems can theoretically support very high data rate over wireless channels without the need of 
additional bandwidth and transmitted power. Thus, the large MIMO systems which comprise of many low power 
antennas could be a future trend for wireless communications that now have very scarce bandwidth. 
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Concatenating the MIMO systems with channel codes is the way to considerably increase the reliability and 
the performance of the MIMO systems. There has been a tremendous effort to develop the coded MIMO systems, 
including code-design, the invention of soft-output detectors, joint detection-decoding techniques and so forth, with 
the ultimate goal to approach the MIMO capacity, (a theoretical minimum SNR at which a reliable communication 
for the given rate is possible), see e.g. (6), Q, (8). However, almost works have studied the coded MIMO systems 
with small number of transmit/receive antennas (e.g. 2 to 8 transmit/receive antennas). So, it is important to note 
that lacking in the literature is a performance study of coded large MIMO system. We also note that the useful 
detection techniques for small MIMO systems, e.g., optimal MAP detection or successive interference cancellation 
techniques, are not applicable to MIMO systems with hundreds of antennas since those methods definitely incur 
ultra high computational complexity at the receiver. 

To the best of our knowledge, only large MIMO systems concatenated with turbo codes have been studied in 1191 , 
ifTOl . ifTTl . [ 12 1 which are the limited number of publications. It is also worth mentioning that the main objective of 
the above cited works is on the development of low complexity large MIMO detections not for an efficient coding 
scheme. So, it is not surprising to see the near optimal detection performance (uncoded performance) but rather 
poor coded performance. For example, turbo coded MIMO system with 200 transmit/receive antennas performs 7 
dB away from the associated capacity which is quite large capacity-gap from the coding theoretic point of view. 
One objective of this paper is to reduce such large capacity-gap by considering the regular NBLDPC codes defined 
over GF(2 8 ) lfT3ll together with the low complexity detections. The reasons why we consider the NBLDPC codes 
will be elaborated in the next section. 

In this paper, we study the application of NBLDPC codes to the MIMO systems with hundreds of antennas that 
utilize low complexity detection algorithms as the MIMO detector. Our main contributions in this paper can be 
summarized as follows : 

1) The NBLDPC coded large MIMO systems with well-known low-compelxity MMSE detector have been 
investigated. By using the soft output generation method proposed in [14], the NBLPDC coded MIMO systems can 
operate nearer to the MIMO capacity comparing with the best known turbo coded systems 0. For example, the 
capacity-gap about 15 dB for turbo coded large MIMO system with 16QAM can be reduced to 6 dB. Moreover, 
under the same soft-output MMSE detection, we have demonstrated via simulation that NBLDPC codes significantly 
outperform both the optimized binary LDPC codes and regular binary LDPC codes over large MIMO systems. 

2) For large MIMO systems, the potential drawback of using MMSE detection is the large size matrix inversion. 
We thus propose the soft output detection based on matched filtering which has lower computational complexity 
(just 0.28% of MMSE detection). Moreover, the similar coded performance, comparing to the case of using MMSE 
detection, in low SNR region or, equivalently, near capacity region can be obtained. 

3) By using the Monte Carlo density evolution lfl5l for analysing the NBLDPC coded large MIMO systems, we 
found that the decoding threshold (minimum SNR at which the decoder can successfully decode the noisy received 
signals) of NBLDPC code is within 1 .6 dB from the MIMO capacity. 

4) Channel estimation is one of the challenges in large MIMO systems. In this study, the sensitivity of the 
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NBLDPC coded large MIMO systems to channel estimation error has been investigated. For both the MMSE 
detection and the proposed detection, large estimation error variance up to 0.2 is still acceptable for NBLDPC 
coded MIMO systems since the loss in coding gain is quite small (only 0.45 dB). 

5) Another one of the main concerns of large MIMO systems is how to place a large number of antennas in 
the limited space. If the antenna spacing is not sufficient, e.g. spacing less than half of wavelength, the MIMO 
systems will suffer from spatial correlated fading which results in the degradation of capacity as well as the coded 
performance. In doubly correlated MIMO channel (correlation at both the transmitter and the receiver), we observed 
that there exist the points in low SNR region at which the effect of spatial correlation is reduced. At this region, 
the NBLDPC coded MIMO system with the proposed detection can be efficiently utilized. 

The rest of this paper is organized as follows. We first describe the NBLDPC coded large MIMO systems in 
Section II. Then the detection problem of large MIMO system is discussed in Section III. In Section IV, NBLDPC 
coded large MIMO systems with MMSE detection are described and investigated. In Section V, we propose a 
new soft output detection for NBLDPC coded large MIMO systems. In Section VI, the thresholds of NBLDPC 
coded large MIMO systems are analysed by Monte Carlo density evolution. The sensitivity of NBLDPC coded 
large MIMO systems is studied in Section VII. In Section VIII, the performance of NBLDPC coded large MIMO 
systems in correlated environment has been investigated. Finally, conclusions are given in Section IX. 

II. System Description 

A. Non-Binary Low-Density Parity-Check Codes 

An NBLDPC code C over Galois field GF(2 m ) |fT3l is defined by the null-space of a sparse P x N parity-check 
matrix A = {a.ij} defined over GF(2 m ), for i = 1, . . . , P and j = 1, . . , ,N 

C = {x e GF(2 m ) N | Ax T = e GF(2" l ) p }, 

where m > 1 and x = (x\, . . . ,xjy) is a codeword. The parameter N is the codeword length in symbol. Assuming 
that A is of full rank, the number of information symbols is K = N — P and the code rate is R = K/N . 

We note that a non-binary symbol which belongs to GF(2 m ) can be represented by the binary sequence of 
length to bits. For each m, we fix a GF(2 m ) with a primitive element a and its primitive polynomial ir. Once a 
primitive element a of GF(2 m ) is fixed, each non-binary symbol is given by an m-bits representation lfl6l p. 110]. 
For example, with a primitive element a e GF(2 3 ) such that ir(a) = a 3 + a + 1 = 0, each symbol is represented 
as = (0,0,0), 1 = (1,0,0), a = (0,1,0), a 2 = (0,0,1), a 3 = (1,1,0), a 4 = (0,1,1), a 5 = (1,1,1) and 
a 6 = (1,0,1). Let L(x) be the binary representation of x E GF(2 m ). For the above example, we can write 
L(x — a 3 ) — (1,1,0). Thus, each coded symbol x t € GF(2 m ),Vi € {1, ...,N} of a non-binary codeword 
represents m bits. We also denote n = mN and k = mK as the codeword length and information length in bit, 
respectively. 

An NBLDPC code is (d v , d c )-regular if the parity-check matrix of the code has constant column weight d v and 
row weight d c . The parity-check matrix A can be represented by a Tanner graph with variable and check nodes ifTTl 
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p. 75]. The belief propagation (BP) algorithm which is commonly used as the NBLDPC decoder |13] exchanges 
the probability vector of length 2 m between variable nodes and check nodes of the Tanner graph at each iteration 
round I. 

In this paper, only is (d v = 2, d c )-regular NBLDPC code defined over GF(2 8 ) considered due to the following 
reasons : 

1) The process to optimize parity-check matrix A is not required since (2, d c ) -regular NBLDPC code defined 
over GF(2 8 ) is empirically known as the best performing code especially for short code length. Moreover, the 
NBLDPC code with d v = 2 can be encoded in linear time ifTHl . 

2) The high decoding complexity of NBLDPC decoder can be compensated since the A of (2, d c ) -regular 
NBLDPC code is very sparse. 

3) For practical code lengths, e.g., a few thousand bits, (2, d c )-regular NBLDPC code defined over GF(2 8 ) 
seems to offer the best performance for small MIMO systems with optimal detection |19|. We therefore expect the 
excellent performance of NBLDPC codes in large MIMO systems. 

4) It is empirically shown that the application of (2, d c ) -regular NBLDPC code defined over GF(2 8 ) to higher 
order modulation is outstanding ll20l p. 32]. 

All good points listed above are very attractive for wireless communications. Thus, it is not overly exaggerated to 
state that we intend to apply the simple, low-complexity, high-performance channel code to large MIMO systems. 

B. System Model 

We adopt the conventional notation to denote the MIMO system with N t transmit antennas and N r receive 
antennas as Nt X N r MIMO system. Let A M be the complex modulation constellation of size M = 2 P where 
p represents bit(s) per modulated symbol. In this study, each antenna uses the same modulation scheme and the 
mapping is a Gray-labelled constellation. We focus our attention to the spatial multiplexing MIMO systems 12T1 . 
E2l that simultaneously send Nt modulated symbols in one time instant (one channel use). At the receiver side, 
some specific techniques will be employed to demultiplex the received data. 

Figure [T] shows the spatial multiplexing MIMO system concatenated with an NBLDPC code. In the case of 
NBLDPC code defined over GF(2 8 ), the system shown in Fig. [I] can be described as follows. At the transmitter 
side, a bit to symbol mapper maps a group of 8 information bits to a non-binary symbol in GF(2 8 ). The stream of 
K symbols in GF(2 8 ) is encoded into a codeword of length N symbols in GF(2 8 ) through an NBLDPC encoder. 
Assuming q > p and q is divisible by p, each coded symbol in GF(2 8 ) is demultiplexed to a group of q = 8/p 
symbols according to the chosen modulation scheme. This means a coded symbol in GF(2 8 ) will be mapped to 
q = 8/p modulated symbols by a constellation mapping. At each time instant, the transmitter simultaneously sends 
N t — K t q modulated symbols in parallel through N t transmit antennas where K t is a number of coded symbols 
per each transmission. 

Let s = [si, S2, ■ ■ ■ , SN t ] T € C Nt be the transmitted vector whose each entry is a complex modulated symbol. 
Each entry Si,Vi € {!,..., N t } taken from A M must satisfy the component-wise power constraint E[||s,;|| 2 ] = 
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Fig. 1. Schematic diagram of the NBLDPC coded MIMO system. For NBLDPC code defined over GF(2 8 ), every 8 bits emitted from source 
will be mapped to a symbol in GF(2 8 ). 



E s /N t where E s is the total transmitted power, || • || 2 denotes the squared Euclidean norm, and E[-] denotes the 
expectation. With this power constraint, a large number of transmit antennas imply low power consumption per 
each transmit antenna. 

Let us consider a 200 x 200 MIMO system with QPSK modulation (M — 4 and p = 2) as an example. 
After encoding, each coded symbol in GF(2 8 ) is mapped to q = 4 modulated symbols. At each time instant, the 
transmitter collects 200 modulated symbols mapping from K t = 50 coded symbols. These 200 modulated symbols 
are simultaneously sent through 200 transmit antennas in parallel. 

Considering a baseband discrete time model for an uncorrected Rayleigh flat fading MIMO channel, the received 



vector y 



[2/1,2/2, ■ ■ ■ ,VN r 



£ C Nr of the spatial multiplexing N t x N r MIMO system is given by |9| 



Hs + n. 



(1) 



The matrix H = [H1H2 . . . Hjy t ] € C rXiv * denotes the channel fading matrix whose each entry h^j is assumed 
to be an i.i.d. complex Gaussian random variable with zero mean and unit variance E[||/ijy|| 2 ] = 1. The vector 



\ni,n 2 , 



■ , n Nr ] T e 



is a noise vector whose entry is an i.i.d. complex Gaussian random variable with 
zero mean and variance er 2 per real component. The MIMO detector performs detection and produces the prior 
probabilities (soft output) for NBLDPC decoder. After all N variable nodes are initialized, the NBLDPC decoder 
performs decoding process and provides the estimated non-binary symbols (hard output). These estimated symbols 
are finally demapped to a sequence of estimated information bits. 
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C. SNR definition and Ergodic Capacity 

In this paper, the channel matrix H is assumed to be known at the receiver and we only focus on the square 
channel matrix, i.e. N t = N r , Since each entry of H has unit variance, the average signal energy per receive 
antenna is E s . We follow the convention that No/2 = <r„ to define the signal to noise ratio. In this setting, the 
average signal to noise ratio (SNR) per receive antenna, denoted by 7, is given by |6] 

E s E s 

7 = tt~ = 7T~2- (2) 

The spectral efficiency (transmitted information rate) of coded MIMO system with spatial multiplexing technique 
is pRN t J5). With perfect H at the receiver side, ergodic MIMO capacity is given by J2) 



where the superscript H denotes the Hermitian transpose operator, det denotes the determinant and I/v r is the 
identity matrix of size N r x N r , Both of spectral efficiency and capacity are measured in bits/sec/Hz (bps/Hz). 

D. Large MIMO Detection 

One of the major problems that prohibits large MIMO systems from the practical realization is the high com- 
putational complexity of the MIMO detector. To clearly show this, let us illustrate by using an example. Suppose 
that the number of transmit antennas are Nt — 100 and the BPSK modulation is employed. A total number of 
possible sequences that can be transmitted from the MIMO transmitter is 2 pNt — 2 100 sequences which are very 
large number. Therefore, performing optimal maximum likelihood detection, i.e., finding the most likely transmitted 
vector given the received signal, on large MIMO systems is infeasible. We also note that the near optimal detection 
algorithms inventing for small MIMO systems, e.g., sphere detection or successive interference cancellation, are 
also too complex for large MIMO systems. Therefore, many works have recently focused on inventing the low- 
complexity detectors to enable the use of large scale multiple antennas 0, (9), iflOl . fl2l . 

In order to operate near capacity region over large MIMO systems, we exactly need the efficient coding and 
detection scheme. Good detection performance (uncoded performance), e.g., near optimal detection performance, 
does not imply good coded performance, e.g., near capacity performance. For example, likelihood ascent search 
(LAS) detection can provide near optimal detection performance but the capacity-gap from using Turbo codes is 
more than 15 dB for 600 x 600 MIMO systems with 16-QAM |9|. This is because, for modern channel codes such 
as Turbo and LDPC codes, good coded performance needs the generation of soft output from the MIMO detector 
as addressed in IflOl . Thus, we believe that the generation of soft output from the MIMO detector is very crucial 
to obtain the good coded performance from using the NBLDPC codes. 

For large MIMO systems, it is well-known that two types of low-complexity linear detections, named as 
minimum mean square error (MMSE) and the matched filter (MF) detections, can provide the near optimal detecting 
performance especially in low SNR region [3|. The results given in Fig. [2] which displays the detecting performance 
of MMSE and MF detections over 800 x 800 MIMO systems with BPSK modulation can be used to confirm the 




(3) 
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above statement. In this figure, the uncoded performance of SISO unfaded AWGN-BPSK system is used to represent 
the approximated lower bound for optimal maximum likelihood (ML) performance (23] . From this figure, the major 
observations could be made : 

• There is no performance difference between MF and MMSE detections when 7 is low, e.g., 7 < —2 dB. 

• The performance of MMSE and MF detections at low SNRs achieves near the performance of SISO unfaded 
AWGN. Thus, both are near optimal in low SNR region. 

• The minimum SNR to achieve capacity below 400 bps/Hz is in the near optimal region of both the MMSE 
and MF detections. So, it is intuitive to think that the MMSE and MF detections with soft output generation can be 
used in conjunction with the NBLDPC codes that have low to medium rates to achieve near capacity performance. 
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Fig. 2. Uncoded performance of MMSE and MF detectors in 800 X 800 MIMO system with BPSK modulation. 



Motivated by the above discussion, we will focus in this work the MMSE and MF detections. An important 
thing that need to be subsequently pointed out is the methodology to obtain the soft output from MIMO detector. 
In the next section, we will present the detection principle of MMSE detection as well as the suitable soft output 
generation that works well with NBLDPC decoder. In Section V, we will develop the novel soft output MF-based 
detection that can provide excellent coded performance with very low computational complexity. 

III. MMSE Detection 

Let us first introduce in this section the MMSE detector and its soft output generation. After that, the coded 
performance will be investigated and discussed. 
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A. Detection Principle 

Following [24], a detection estimate Sfe of the transmitted symbol on fc-th antenna by multiplying with the 
MMSE weight vector W fe 

h=W%y k , (4) 
where the MMSE weight matrix Wfe is of the form 

W fe = ( -J^I Nr + HH H ) H A . (5) 



K E s /N t ' N ' 

This MMSE weight vector Wfe is chosen so as to minimize the mean square error between the transmitted symbol 
Sfe and ife. 

B. Soft Output Generation 

By using the method for obtaining soft output from MMSE detector given in IPT41 . the estimation Sfe can be 
approximated as the output of an equivalent additive white Gaussian noise (AWGN) channel 

h = MfcSfe + zh, (6) 

where p,f. = WjjHfc and is a zero-mean complex Gaussian variable with variance 4 = jf(p,k — A*fe)- Based on 
this approximation, the probability of Sfe conditioned on s E A M is as follows 

Pr(J fe | s) ~ -^exp ( -- j || s fe -fi k s || 2 ) , (7) 
*4 \ 4 J 

where the parameter n is the normalized constant such that J2 s ea m ^ r I s ) = P° r 171 7^ P which is the 
condition to be considered, the probability Pr (Jfe | s) cannot feed directly to NBLDPC decoder. 

By deploying the iterative belief propagation (BP) algorithm as an NBLDPC decoder J25), the procedure to 
form the input of NBLDPC decoder which is defined over GF(2 8 ) from the known probabilities calculating from 
is as follows. Let p^\x) denote the probability that v-th coded symbol is likely to be x E GF(2 8 ) where 
v = 1,2, ... ,N and the index in superscript represents the iteration round of BP algorithm. We assume that fc-th 
transmit antenna to (fc + q — l)-th transmit antenna are used to send the v-th coded symbol (any coded symbol 
is mapped to q modulated symbols). For each x E GF(2 8 ), we need to collect q modulated symbols to represent 
qp = 8 bits. Thus, the generation of soft output from MMSE detector for the symbol x E GF(2 8 ) of v-th variable 
node can be expressed as 

9-1 

pM(x) = l[Px(s k+ i\s k+i ), (8) 

i=0 

where Pr (sfe+i | s k+i ) is the likelihood of the estimation Sk+i conditioned on s k+t E A M by which the eight 
bits ordered sequence of (s k , s k+1 , ■ ■ ■ s k+q ~ 1 ) must be equal to the binary representation L(x). An example of 
calculating pffl (x) is illustrated in Fig. [i] 

According to 2 s (q — 1) real multiplications for each coded symbols are needed to calculate p^\x) for all 
x E GF(2 8 ). We note that the computational complexity of generating soft output for NBLDPC decoder is rather 
low when comparing to the computation of MMSE matrix. 
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Fig. 3. The example of calculating soft output for symbol x with L{x) = [00 01 11 10]. In this example, we assume that the complex 
modulation constellation of size M = 4, e.g., QPSK, is employed and the elements in this constellation are a, b, c, d with their gray-labelled 
representations 00, 01, 11, 10, respectively. 

C. Performance 

In this part, the performance of NBLDPC coded MIMO systems with soft output MMSE detector is presented. 
The maximum iteration of NBLDPC decoder £ max is set to 200 for all simulation results |26|. 

The BER performance of NBLDPC coded 200 x 200 MIMO systems with BPSK modulation (M = 2) is shown 
in Fig. [4] For 200 x 200 MIMO system, the best performing scheme which can be founded in the literature is 
the R = 1/2 turbo coded spatial multiplexing MIMO systems with the MMSE-LAS detector 0. The MMSE- 
LAS detection algorithm uses the MMSE detection to initialize the algorithm. Therefore, the overall computational 
complexity of MMSE-LAS detector is greater than that of MMSE detector. The BER performance of this turbo 
coded system is away from MIMO capacity by 7.5 dB. It is clearly seen from the figure that R = 1/2 NBLDPC 
coded system significantly outperforms turbo coded system by about 4 dB. Both R = 1/3 and R— 1/2 NBLDPC 
coded systems perform within just 3.5 dB from the MIMO capacity. 

Figure [5] presents the simulation results for larger dimension, 600 x 600 MIMO system. Comparing to 200 x 200 
MIMO system, it can be observed that the capacity-gaps remain the same (3.5 dB) for both R = 1/3 and R = 1/2 
NBLDPC coded systems. However, the capacity-gap of turbo coded system with MMSE-LAS detector in |9| is 
larger (about 9.4 dB) when the dimension of channel increases. It is obviously seen that R — 1/3 turbo coded 
system underperforms R = 1/3 NBLDPC coded system by more than 6 dB. More interestingly, R = 1/2 NBLDPC 
coded system absolutely outperforms R = 1/3 turbo coded system by about 3.4 dB. Therefore, the NBLDPC coded 
system outperforms the turbo coded system both in terms of performance and spectral efficiency. 

It is well known that the spectral efficiency can be enhanced without additional bandwidth by adopting higher 
order modulations such as M-QAM or M -PSK. Thus, it is interesting to investigate the performance of coded large 
MIMO system with higher order modulation. The coded performance of 600 x 600 MIMO system with 16-QAM 
is shown in Fig. [6] From this figure, the following observations can be listed as follows : 

• R = 1/3 and R = 1/2 NBLDPC coded systems with n = 2400 bits operate within 6 dB and 8 dB, respectively, 
from MIMO capacity. 
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Fig. 4. Bit error rate performances of coded 200 X 200 MIMO systems with BPSK modulation. The spectral efficiencies of 66.6 and 100 
bps/Hz are obtained from MIMO systems with channel codes of R = 1/3 and R = 1/2 respectively. 



• R = 1/3 and R = 1/2 turbo coded systems |9j operate very far (more than 15 dB) from MIMO capacity. It is 
clearly seen that the NBLDPC coded systems indeed beat turbo coded systems. 

• By increasing the length of code from n = 2400 bits to n — 28800 bits, the coding gain about 0.7 dB can be 
obtained from R — 1/3 NBLDPC coded system. 

Based on the result, we can conclude that the good performance over large MIMO channels by using NBLDPC 
codes and 16-QAM are possible. 

One may think that the previous results are not sufficient to show the benefit of using NBLDPC code in large 
MIMO systems. Moreover, we conduct the comparison among different types of LDPC codes : regular NBLDPC 
code over GF(2 8 ), optimized binary LDPC code, and regular binary LDPC code. By fixing the code rate R = 1/3, 
the regular binary LDPC code with column weight 4 and row weight 6 and the optimized binary LDPC code from 
Table. Ill in lf27l are chosen as the competitors. The maximum degree of variable node of the selected optimized 
binary LDPC code is 16. We choose this optimized LDPC code because it can asymptotically perform very close 
to the capacity of single input single output Rayleigh fading channel (within 0.19 dB). By using the same 16-QAM 
and soft output MMSE detection, the performances of 600 x 600 MIMO systems coded with different kinds of 
LDPC codes are shown in Fig. [6] From this figure, the interesting observation can be listed as follows : 

• R = 1/3 NBLDPC coded systems outperform both the optimized and regular binary LDPC (BLDPC) codes. The 
coding gains obtained over optimized and regular binary LDPC codes are about 0.8 and 2 dB respectively. 

• Another advantage of using NBLDPC code which can be seen from this figure is the excellent frame error rate 



11 



10 



c 

h 
m 



10 



— min. y for capacity = 200 bps/Hz 
' — min. y for capacity = 300 bps/Hz 
-4r— R-l/3 Turbo-MMSE-LAS [9] 
-a— R-l/3 NBLDPC-MMSE (n=2400 bits) 
e - R-l/2 NBLDPC-MMSE (n=2400 bits) 




-5 -4-3-2-1 1 2 3 

Y : Average SNR per Receive Antenna (dB) 



Fig. 5. Bit error rate performances of coded 600 X 600 MIMO systems with BPSK modulation. The spectral efficiencies of 200 and 300 
bps/Hz are obtained from MIMO systems with channel codes of R = 1/3 and Ft = 1/2 respectively. 



(FER) which is an advantage that can be expected from LDPC codes of column weight two |28|. Comparing with 
optimized BLDPC code, the FER of NBLDPC code is much better. Although the BER curve of optimized BLDPC 
code is good but its corresponding FER curve is quite bad. 

We have shown in this section that NBLDPC coded large MIMO systems with soft output MMSE detection 
can be employed to reduce the remaining capacity-gap of the best known turbo coded systems. Despite the low 
complexity of MMSE detection, one may think that finding inverse of large and dense channel matrix which is 
the main task of MMSE detection is rather cumbersome. Inspired by the very low-complexity and near optimal 
detecting performance of MF detector, the soft output MF-based detector will be developed in the next section. 



IV. The Proposed MF-Based Detection 

A. Detection Principle 

A baseband discrete time channel model given in ([1} can be equivalently rewritten as 

y = Hisi + H 2 s 2 H h HN t SN t + n. 



(9) 



To estimate §k of the transmitted symbol on the k-th antenna by the concept of matched filtering, the interference 
streams in |9) are treated as noise, 



y = H fe s fc + ^2 H * s « + n ■ 



(10) 
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Fig. 6. Bit error rate performances of coded 600 X 600 MIMO systems with 16-QAM modulation. The spectral efficiencies of 800 and 1200 
bps/Hz are obtained from MIMO system with channel codes of R = 1/3 and R = 1/2 respectively. 



Then, the MIMO detector based on MF technique estimates §k of the transmitted symbol on the fc-th antenna by 

jjH 

multiplying the received vector y with a vector Wfe 



H H t 



Sfe = W fe y, 

N t (11) 
= s k + Wfe £ H lSi + Wfen. 

i—l^i^k 

In contrast of the conventional MF detection, we make a little modification by introducing the term ttH 1 tt into 

H fc H fc 

weight vector. After slicing Sfe, we can obtain the estimated symbol transmitted from the feth antenna. It is known 
that the detection performance (uncoded performance) of match filtering is very poor since there is no cancellation 
of the interference terms (the second term in R.H.S of ([TT|)). However, matched filtering is near optimum if n is 
dominant Q, 11291 , i.e., at low SNRs, as shown in Fig. [2] 

To reduce the computational complexity of detection, we observe that E [H^Hfe] w N r . Thus, in the case of 
hundreds of antennas, the computation of Wfe can be simplified as follows 

Wfe =H»{1/N r ). (12) 

With this simplification, the matrix multiplication H^Hfe whose an complexity is 0(N r ) ll30l can be reduced to just 
a constant term. Figure [8] shows that there is no performance difference between MF detection and its corresponding 
simplified version. 
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Fig. 7. Performance comparison between R = 1/3 NBLPDC code and R = 1/3 BLDPC codes on 600 X 600 MIMO systems with 16-QAM 
modulation at n = 2400 bits. The solid curves represent the BER performance while the dashed curves are the corresponding FER performance. 



B. Soft Output Generation 

In this part, a novel soft output generation from MF-based detection for NBLDPC decoder which has not 
yet reported elsewhere is presented. Assuming all streams and noise are statistical independence and Gaussian 
distributed, the proposed method for generating soft output from MF-based detection will be described as follows. 
Following the definition of the post detection signal to interference plus noise ratio (SINR) given in OT1 p. 358], 
the SINR of the proposed detection for k-th stream, denoted by 8k, can be expressed as 

5, = E IIMI 2 1 

E[ J2 ||W fc H iSi || 2 ] + E[||W fc n|| 2 ] 

= ^Bi (13) 

Es/Nt E l|W fc H s || 2 + (2<r 2 ||W fc || 2 ) 
_ E 3 /N t 



The denominator, denoted by Afc, given in ( 13 i can be approximated by a Gaussian random variable. Based on 
this approximation, we set u\ = A&/2 and the soft output which exactly is the likelihood of s& conditioned on 
s € A A/ is as follows 

Pls - ]s}= 7m^(-^ h -*)- <14) 

To justify that distributes like Gaussian, we estimate the probability density function (pdf) of random variable 
Afe by mean of Kernel density estimation |32|. By using 100,000 channel realizations, it is obviously seen from 
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Fig. 8. Uncoded performance of MMSE, ZF, MF, and simplified MF detectors in 200 X 200 MIMO system with BPSK modulation. The 
uncoded performance of SISO unfaded AWGN-BPSK system is again used to represent the approximated lower bound for ML performance 
1231 . The performance of zero forcing (ZF) detection which is a well-known low-complexity detector is also shown. The ZF detection performs 
very poor in low SNRs so we do not pay the attention to this type of detector. 



Fig. [9] that the pdf of A/, closely agrees with the shape of Gaussian distribution. Moreover, we also compare 
the distribution of A& and Gaussian distribution by means of Kolmogorov-Smirnov test |[33l pp. 623-626]. It is 
found that both have the same distribution at the 0.1% significance level. Therefore, we conclude that the Gaussian 
approximation of A^ described above is reasonable. By using (|8j, the input for NBLDPC decoder can be now 
formalized from ( fl4) , 

To reduce the computational complexity of soft output generation, we first observe that E [||W fe || 2 ] « 1/N r . 
This is because E [||H^!|| 2 ] w N r and E [||H^Hfc|| 2 ] sa iV 2 . Then, we assume that the power of noise term is 
dominant, i.e., 2tr 2 3> E s . With this assumption, we have the following expression 

E s /N t l|W fc H J || 2 «2cT2||W fc || 2 . 

i— l : i^k 

N t 

The above expression is valid when 7 is low, e.g., 7 = -5 dB, since l/N t Yl ||WfcH,;|| 2 < ||Wfc|| 2 . Therefore, 

i— l,2^fc 

the denominator term A^. in (13i can be simplified as follows, 



A fc =E S /N t £ HWfcHJ 2 + (2a*\\W k 



<k\ 

r 2 



(15) 



2al/N r . 



With this simplification, the computation of variance a\ for soft output generation defined in (I61 can be reduced 
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Fig. 9. The pdf estimation of A fc from 200 X 200 MIMO channel with BPSK modulation at SNR per receive antenna 7 = -2 dB. 



to a constant term aj, — Afe/2 = a^/N r . The variance a\ for soft output generation is now independent of k 
and can be easily pre-computed. We note again that the method for generating soft output from MF-based detector 
presented in this part is completely different from the previously found techniques for obtaining soft output from 
MF detection l34l. l35l. 



C. Performance 



Figure 10 shows the performance comparison between the NBLDPC coded systems with the proposed soft output 
MF-based detector and the NBLDPC coded systems with soft output MMSE detector described in the previous 
section when the number of transmit antennas is N t — 200. The rates of NBLDPC codes are ranging from 1/12 
to 1/3 and the information length is 800 bits (K = 100 symbols). As mentioned in Section III, we select the low 
rate NBLDPC codes since we expect that the operating region, e.g. BER of 10 -4 , will be occurred in low SNRs at 
which the proposed detection is near optimal. For R = 1/3, the performance of coded system with MMSE detector 
is better than that of the proposed system by about 0.8 dB. This is because the uncoded performance of MF-based 
detector is worse than that of MMSE detector by roughly 2 dB (not show here). However, as we expected, the 
performance difference between two coded systems (one with MMSE detector and another one with the proposed 
detector) is vanished for R = 1/6,1/9 and 1/12 since, in low SNRs, both detectors provides almost the same 
detecting performance. Thus, the soft output generation from MF detector presented in this paper works well with 
NBLDPC decoder. Note that the codes with R < 1/3 are constructed from NBLDPC of rate 1/3 by using the 
instruction given in |36|. 
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Fig. 10. BER of the NBLDPC coded 200 X 200 MIMO systems with BPSK modulation. "Propose" means the NBLDPC coded system with 
the proposed soft output MF-based detection and "MMSE" represents the non-binary LDPC coded system with soft-output MMSE detection. 



We have shown in Fig. 10 that the application of the proposed soft output MF-based detector in NBLDPC coded 
MIMO system with binary modulation is excellent if the rate of NBLDPC code is low, e.g., R < 1/3. The next step 
is to illustrate the coded performance with higher order modulation. Figure [TT] shows the performance of NBLDPC 
codes of R = 1/9, 1/12 when the 16-QAM modulation is employed. Unlike the case of BPSK modulation, we 
can see the error floor at BER of 10 -4 . The error floor appears because the proposed detection has no interference 
cancellation which is more severe to the case of higher order modulation (more points in modulation constellation 
comparing with the binary case). However, we observed that the number of bit errors in each frame error is small. 
For example, 100 frame errors at 7 — 2 dB for NBLDPC code of R = 1/9 equal to 1752 bits error. The average 
bit errors per one frame error is just 18 bits which is rather small. We intend to use the joint detection-decoding to 
overcome the error floor but this topic is left as the future work. 

In this part, we have demonstrated the application of the proposed soft output MF-based detector in NBLPDC 
coded large MIMO systems. The proposed detector can be efficiently used with both binary modulation but we 
mention that this kind of detector is useful in the limited operation region, i.e., the region that MF detector is 
near optimal. In that region, we need to use the low rate channel code and the low rate NBLDPC codes are the 



promising candidate as we shown in Fig. 10 Please note that the NBLDPC codes of R < 1/3 used in this work 
can be decoded at almost the same computational complexity [36 1 by using the Tanner graph of NBLDPC code of 
R = 1/3 which is another advantage of using NBLDPC codes. 
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Fig. 11. BER of the NBLDPC coded 200 X 200 MIMO systems with the proposed soft output MF-based detector and 16-QAM modulation. 
The information size k = 800 bits and the code length is 800/i? bits. 



D. Complexity Analysis 

Following 071 . If30l p. 18], the complexity of complex matrix calculation is expressed in terms of complex 
floating point operations (flops). The number of flops for some mathematical operations is given in Table [I] We 
note that the detection for coded MIMO systems comprises of two steps : 1) detection and 2) soft output generation. 



TABLE I 

The number of flops for some operations f37l and |[38] p. 193] for exponential calculation where a, b e K, a, b e C, and 

c,d 6 C N - X1 



Operation (name) 


Operation (math.) 


Number of flops 


Real multiplication 


ab 


1 


Complex multiplication 


ab 


3 


Real addition 


a + b 


1 


Complex addition 


a + b 


1 


Inner product 


c H d 


AN r - 1 


Scalar-vector multiplication 


a • a 


N r 


Exponential calculation 


exp(a) 


50 
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1 ) Proposed Soft Output MF-based Detector K39\l : The estimation of §k as defined in ( 1 1 1 needs two steps : 1) 



Computing the weight matrix according to ( [12} which requires N r flops. 2) Multiplying with received 
vector y which requires 4N r — 1 flops. Since we have N r receive antennas, the total flops used for detection are 
5N£ — N r . Thus, the per-bit complexity is 0{N r ). To obtain the soft output defined in (14 1, we need one flop 



for computing squared euclidean norm, one flop for subtraction, 3 flops for multiplication of real constant, Each 
transmitted symbol has M levels to be calculated and we have N t = N r transmitted symbols. Therefore, in total, 
55MN r flops are required for generating soft outputs. In summary, the overall per-bit complexity for both detection 
and generating soft output is 0(N r ). 

2) Soft Output MMSE Detector: In the case of MMSE detector, it has been shown that the complexity for 
detection is ION, 3 + b.bNf. + l.5N r flops J33. To generate the soft outputs as defined in fi k = Wj^, fc H fc 
requires AN r — 1 flops and = jf^i^k — Mfc) requires 3 flops. Again, each transmitted symbol have M level to 
be calculated and we have Nt — N r transmitted symbols. So, AMN^ + 58MN r flops are required for generating 
soft outputs from all estimated symbols. 

Table \U\ summarises the computational complexity of both the proposed detector and the MMSE detector in terms 
of flops. In the case of BPSK modulation (M = 2), the computational complexity versus the number of receive 



antennas is plotted in Fig. 12 It is seen from the figure that the computational complexity of the proposed detection 
is much lower than that of MMSE detection For instance, 80, 563, 500 flops are required for the MMSE detector 
while we need only 221,800 flops for the proposed MF detector (just 0.28 % to the MMSE detection) Thus, by 
using the proposed detector, the computational complexity of large MIMO detection is significantly reduced. 

TABLE II 

The number of flops for MIMO detection. 



Type \ Operation 


Detection 


Soft Output 


Proposed detector 


5N? - N r 


55MAT r 


MMSE detector 


ION? + 5.5JV-? + 1.5AV 


4MJV"2 + 58MN r 



V. Threshold of Non-Binary LDPC Codes in Large MIMO Systems 

In the previous two sections, the performances of NBLDPC coded large MIMO systems with both the MMSE 
detection and the proposed MF-based detection have been investigated by using the practical code lengths, e.g., a 
few thousand bits. We found that NBLDPC coded systems have the capacity-gaps about 3.5 dB and 6 dB, which 
can be considered as the best performing for BPSK modulation and 16-QAM, respectively. The fact from coding 
theory tells us that the capacity-gap can be reduced by using the larger code length as demonstrated in Fig. [6] So, 
one may think about the performance of NBLDPC codes of large code lengths. We now show in this section the 
capacity-gap of NBLDPC coded large MIMO systems which have very large code length. 
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Fig. 12. The computational complexity of the MMSE detection and the proposed detection versus the number of receive antennas. 



In this section, we analyse the performance of NBLDPC coded large MIMO systems in the limit of infinite code 
length by using the Monte Carlo density evolution proposed in fl5l . The Monte Carlo density evolution is a tool 
for calculating the minimum SNR (known as decoding threshold) at which the NBLDPC codes can reliably decode 
the noisy received symbols. The calculating steps of Monte Carlo density evolution can be summarized here. 

1) Input the size of variable node L, number of iteration £ max , and the initial SNR 7^ in dB. The size of variable 
node L must be large enough to represent the infinite code length, e.g., 10 5 nodes. The maximum iteration £ max is 
also large, e.g., 2000 iterations. The initial SNR 7^ should start from the value that NBLDPC decoder can reliably 
decode., e.g., the SNR that can provide BER of 10~ 4 in the case of practical code length. Let £ be the index of 
iteration. 

2) By assuming the zero codeword to represent the NBLDPC codeword, we first produce L noisy variable nodes 
with associated likelihoods p^\xj) from the equivalent channel as shown in Fig. 

Ij eGF(2 m ) ! j = i,..,r. 

3) Increment £ — £ + 1. Creating the new set of variable nodes of size L with the associated likelihood pi\x) by 
using the concept of BP decoder [13]. The operation in this step, which is illustrated in Fig{l4] is like running BP 
algorithm on cycle-free Tanner graph. The full detail of this step can be found in [15 p. 22]. 

4) Replacing the previous set of variable nodes of size L with the new set of variable node of the same size. 
Calculating the average Shannon entropy as shown below 

H = -\ (^) lo g 2m (pi e) (xjj) , (16) 

v=l j=l 



13 



where v = 1,...,L and 
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where, in this case, the value of H can be used to represent the ambiguity of the underlying variable nodes. Thus, 
the value of H = implies a complete removal of the ambiguity. In this paper, we choose -ff s top = 1CP 6 as the 
stopping criterion. If H > H stop and I < l max , repeat step 3 and step 4 in iterative manner until I = £ max . If 
H < H sto p, go to step 5. If H > H stop and I = £ max , go to step 6. 

5) If H < H stop , this means the NBLDPC decoder can reliably decode the noisy variable nodes at 7W. Decrease 
the 7^ by small step size, e.g., 7^ +1 ) = jW — 0.05 dB and then repeat step 2 to step 4. 

6) If H > -ff s top and I = £ max , this means the NBLDPC decoder cannot reliably decode the noisy variable nodes 
at 7W. Declare the previous as the decoding threshold. 
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Fig. 13. An equivalent channel of NBLDPC coded MIMO systems. 




Fig. 14. The generation of new set (ensemble) of variable nodes in the Monte Carlo analysis originally proposed by M. Davey 1151 . By 
iterating the computation to certain iteration (quite large), this is equivalent to perform decoding on cycle free Tanner graph with infinite size 
of variable nodes. 



Figure 



15 



shows the decoding threshold of NBLDPC coded 200 x 200 MIMO system with the proposed MF-based 
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detection and BPSK modulation. We also plotted the MIMO capacity of 200 x 200 MIMO system and the SNR 
points required to reach the BER of 10~ 4 of the NBLDPC coded system with short information length k = 800/i? 
bits. As you can see from the figure, the capacity-gap of NBLDPC codes does not go to small value even when 
the code length tends to infinity. The capacity -gap of NBLDPC codes to the MIMO capacity is about 1.6 dB for 
all rates. Therefore, we cannot achieve the MIMO capacity by simple concatenation the large (spatial multiplexing) 
MIMO systems with regular NBLDPC codes over GF(2 8 ). The methodologies to reduce the remaining capacity-gap 
are challenging and will be explored in the future work. 

We also confirm the correctness of the decoding threshold obtained from Monte Carlo analysis by conducting 



the simulation of the code with very long code length. It can be observed from Fig. 16 that the frame error rate 
converges to the decoding threshold as the code length increases. 




y : Average SNR per Receive Antenna (dB) 



Fig. 15. The decoding threshold of NBLDPC coded 200 X 200 MIMO system with BPSK modulation. The corresponding MIMO capacity 
and the performance at code length n = 800/ R bits are also plotted. 



We note that the same Monte Carlo density evolution cannot be used to compute the decoding threshold of 
NBLDPC coded MIMO systems with M-QAM modulation where M > 2. This is because the M-QAM modulation 
does not has the property of rotational symmetry so the zero codeword cannot be used to represent the LDPC 
codeword. The modified Monte Carlo density evolution presented in |40| can be used to compute the decoding 
threshold of NBLDPC coded MIMO systems with M-QAM. However, this method is computational intensive so 
we leave the decoding threshold of NBLDPC coded MIMO systems with M -QAM as the future work. 
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Fig. 16. Frame error rate of NBLDPC codes of R = 1/6 over the 200 MIMO channel with BPSK modulation. 



VI. Effect of Channel Estimation Errors 

Beside high detection complexity, channel estimation is also a major bottleneck in large MIMO systems 1411 . 
In this section, the effect of imperfect channel estimation on the NBLDPC coded large MIMO systems will be 
investigated. After performance evaluating, we expect to answer the readers about the sensitivity of NBLDPC 
coded large MIMO systems to the accuracy of channel estimation. By using the model adopted in [9|, the estimated 
channel matrix to be used at the receiver is given by 

H = H + AH, (17) 

where AH is the estimation error matrix in which the entries are assume to be the complex Gaussian with zero 
mean and variance ai. 



Figure [TTj shows the BER of the NBLDPC coded 200 x 200 MIMO systems with BPSK modulation for different 
values of a 2 e — 0, 0.1, 0.2. For both the MMSE and the proposed detections, the coded performance degrades when 
(Tg increases as expected, e.g., worse channel estimation quality, worse coded performance. However, the figure 
clearly shows that the performance degradation is negligible even at large error variance <j\ = 0.1 (the error variance 
at this level is considered to be large ll42l ). Note that the performance with a 2 < 0.1 is identical to that of a\ = 0.0 
so we do not included in the figure. For larger error variance a 2 — 0.2, the performance degradation is only 0.4 
dB which is still small. The same property can also be seen although the channel dimension is increased from 200 



to 600 as shown in Fig. 18 



The channel estimation algorithm that can render the error variance to be 0.2 is enough to provide not too large 
performance loss. The reasons why the effect of channel estimation errors to NBLPDC coded large MIMO systems 
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is quite small are : 1) the effect of channel estimation errors for uncoded performance in near capacity region 
is small, 2) the powerfulness of NBLDPC codes, 3) if one look at the simplification of the proposed MF-based 
detection, one may think that this detection scheme does not require the perfectness of channel information at 
receiver side. Based on the results presented in this section, it is sufficient to conclude that the NBLDPC coded 
large MIMO systems with MMSE and proposed detections are very robust to channel estimation errors. 




-4 -3.5 -3 -2.5 -2 -1.5 

y : Average SNR per Receive Antenna (dB) 



Fig. 17. BER of the NBLDPC coded 200 X 200 MIMO systems with BPSK modulation for different values of cr^ by which crj? = means 
perfect channel estimation. In this figure, the NBLDPC code of R = 1/3 and n = 2400 bits is employed "Propose" means the non-binary LDPC 
coded system with the proposed soft output MF-based detection and "MMSE" represents the non-binary LDPC coded system with soft-output 
MMSE detection. 



VII. Effect of Spatial Correlation 

All the previous results of NBLDPC coded large MIMO systems have been evaluated by assuming i.i.d. fiat 
fading channel, i.e., each entry of H is modelled as identically and independent complex Gaussian random variable. 
However, the realistic MIMO channels have spatial correlation, i.e., the entries of H (the propagation coefficients) 
are correlated, which has the adverse effect on the MIMO capacity as well as the error rate. It is known that low 
angle spread, insufficient antenna spacing and lack of rich scattering may cause the spatial correlation between 
antennas ||431 . For large MIMO systems, the insufficient antenna spacing is an important factor to cause the spatial 
correlated MIMO channel since a large number of antennas, e.g., 200 transmit/receive antennas, need to be placed 
in a limited space. So, the study of spatial correlation in NBLDPC coded large MIMO systems should be addressed. 

Since, in this work, both the transmitter and the receiver have large number of antennas, we need to consider 
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Fig. 18. BER of the NBLDPC coded 600 X 600 MIMO systems with BPSK modulation for different values of o*. In this figure, the NBLDPC 
code of R = 1/6 and n = 4800 bits is employed "Propose" means the NBLDPC coded system with the proposed soft output MF-based 
detection. 



the doubly correlated MIMO channels models which can be written as [44 1 



(18) 



where Hi.i.d. is the channel fading matrix of size 7V r x N t whose entry is assumed to be i.i.d. complex Gaussian 
random variable with zero mean and unit variance, R r is the received correlation matrix of size N r x N r , and R 4 
is the transmitted correlation matrix of size N t x N t . 

In fact, there are many ways to generate the correlation matrices R r and R t . To study the effect of correlation 
on the MIMO systems, we first consider an simple and well-known exponential correlation model as shown below 



Rr {Pr) 

Rt (p t ) 



Ji-i\ 

Pr 
Pt 



i,j=l,...,N r 



Pr G [0,1), 

, Pt G [0, 1), 



(19) 



i,j=l,...,N t 

where p t and p t denotes the transmitted and received correlation parameter, respectively. This model is reasonable 
for equally-spaced linear antenna array [44|. Furthermore, this model is physically reasonable in the sense that the 
correlation decreases with the increasing of distance between transmit/receive antennas |451 . 

Figure 19 shows the effect of correlation on MIMO capacity when N t = N r — 600 and the modulation is BPSK. 
We denote parameter p = p t = p r as the correlation parameter which is the same for both the transmitter and the 
receiver. The parameter p — means uncorrected system whereas p = 1 means fully correlated system. It is not 
surprising to see that the capacity of MIMO systems decreases with the increasing of correlation parameters. This 
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is because the spatial correlation has the adverse effect on rank structure of H which directly relates to the MIMO 
capacity. Increasing the number of receive antenna is found to be the solution to keep the MIMO capacity as well 
as the coded performance at the same level [10|. 

Fortunately, one can expect small performance degradation in low SNR region since the effect of correlation in 
low SNR region is quite small. For example, the capacity at 7 — — 11 dB for uncorrected case is about 64 bps/Hz. 
We observe that the loss in capacities at 7 = —11 dB for p = 0.3, 0.4, 0.5 are just 0.85, 1.5, 2.7 bits, respectively, 
which are rather small if one look at the loss in high SNR region. At low SNR region, we know that the proposed 
MF-based detection which has very low computational complexity can be efficiently utilized. To illustrate this, we 
show the performance of NBLDPC coded 600 x 600 MIMO systems at spectral efficiency 66.67 bps/Hz in Fig. 
20 Comparing to uncorrelated case p = 0, the performance loss from p — 0.3, 0.4, 0.5 are approximately 0.4,1,1.5 
dB, respectively. If the bandwidth is 15 MHz (not too large bandwidth requirement), the NBLDPC coded large 
MIMO systems that suffer from the spatial correlation could be the promising candidate to provide the the reliable 
communication (BER of 10~ 4 ) with data rate 1 Gbps (very high) and the received power around -8 dB (very low). 
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Fig. 19. Capacity of 600 X 600 MIMO system by using an exponential correlation model with different correlation parameters p ■■ 



Pt- 



Fig. 21 shows the MIMO capacity under correlated environment modelled by the instruction given in II46II . This 
model allows us to study the effect of antenna spacing which is of interest in large MIMO systems. The MIMO 



capacity decreases if the antenna spacing decreases as seen from Fig. 21 Like the exponential correlation model 



that we have previously discussed, there still exists the points in low SNR region, at which the effect of spatial 
correlation is rather small. So, at that region, we can utilize the NBLPDC coded system with the proposed ultra 
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Fig. 20. Performance of R = 1/9 NBLDPC coded 600 X 600 MIMO systems in the presence of correlation (exponential correlation model). 
BPSK and the proposed detection are utilized as modulation scheme and detector, respectively. The spectral efficiency is about 66.67 bps/Hz 
and the information length is k = 800 bits. 

low-complexity MF-based detection to construct the excellent MIMO communication although the antenna spacing 
is less than half of wavelength. 

VIII. Conclusions 

In this work, we have demonstrated the excellent performance of (2, d c )- regular NBLDPC codes over GF(2 8 ) 
in the MIMO systems with hundreds of antennas. The advantage of using NBLPDC codes in these large MIMO 
systems can be listed here : 1) very good error rate by using the low-complexity detections, 2) robust to channel 
estimation errors, 3) application in correlated environment is excellent. We also proposed in this work the ultra 
low-complexity soft output MF-based detection which can be efficiently utilized for NBLDPC coded systems in 
near capacity region. Supremely, we believe that the results presented in this work can be used as the benchmark 
for the further study on coded large MIMO systems. 
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